Modeling the light curve of the transient SCP06F6 



Emmanouil Chatzopoulos^, J. Craig Wheeler ^, J. Vinko^'^ 

manolisOastro . as .utexas . edu, wheelQastro . as .utexas . edu, 
vinkoQastro . as . utexas . edu 

ABSTRACT 

We consider simple models based on core collapse or pair-formation super- 
novae to account for the light curve of the transient SCP06F6. A radioactive 
decay diffusion model provides estimates of the mass of the required radioactive 
nickel and the ejecta as functions of the unknown redshift. An opacity change 
such as by dust formation or a recombination front may account for the rapid 
decline from maximum. Within this class of model, the redshift must be less 
than 2 ~ 1 or the nickel mass would exceed the total mass of the ejecta, the 
radiated energy would exceed the kinetic energy and kinematic and photometric 
estimates of the radius would disagree. We particularly investigate two specific 
redshifts: z = 0.143, for which Gaensicke et al. (2008) have proposed that the 
unidentified broad absorption features in the spectrum of SCP06F6 are C2 Swan 
bands, and z = 0.57 based on a crude agreement with the Ca H&K and UV iron- 
peak absorption features that are characteristic of supernovae of various types. 
For the lower redshift, we obtain a nickel mass of 0.3 Mq and an ejected envelope 
mass of ~ 38 M©, while for the latter case we find 4.8 Mq and 20 Mq, respec- 
tively, for fiducial parameters. The kinetic energy of the ejecta, while dependent 
on uncertain parameters, is generally large, ~ 10^^ erg, throughout this range 
of redshift. The ejected masses and kinetic energies are smaller for a 
more tightly constrained model invoking envelope recombination. We 
also discuss the possibilities of circumstellar matter (CSM) shell diffusion and 
shock interaction models. In general, optically-thick CSM diffusion models can 
fit the data with the underlying energy coming from an energetic buried super- 
nova. Models in which the CSM is of lower density so that the shock energy is 
both rapidly thermalized and radiated tend not to be self-consistent. We suggest 
that a model of SCP06F6 worth futher exploration is one in which the redshift is 
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~ 0.57, the spectral features are Ca and iron peak elements, and the light curve 
is powered by the diffusive release of a substantial amount of energy from nickel 
decay or from an energetic supernova buried in the ejecta of an LBV-like event. 

Subject headings: circumstellar matter - stars: evolution - supernovae: general 
- supernovae: individual: SN SCP06F6 - hydrodynamics 

1. INTRODUCTION 

Recently discovered very luminous supernovae, such as SN 2006gy (Quimby 2006; Smith 
et al. 2007), SN 2005ap (Quimby et al. 2007a), SN 2006tf (Quimby Castro & Mondol 2007; 
Quimby ct al. 2007b; Smith et al. 2008) and SN 2008es (Yuan ct al. 2008; Gezari et 
al. 2008; Miller et al. 2008) introduce evidence for new supernovae phenomena. The 
discovery of SN 2006gy triggered a rich discussion of the nature of this event and a number 
of models were proposed to interpret the observed light curve (Ofek et al. 2007; Smith et 
al. 2007; Smith & McCray 2007; Woosley, Blinnikov & Heger 2007; Angoletto et al. 2008). 
SN 2006gy is classified as a Type Iln event that peaked at V~22 70 days after the explosion 
and radiated away more than lO^^ergs of energy. Radioactive ^^Co decay fits of the light 
curve of SN 2006gy imply extraordinary amounts of initial nickel mass of the order of 22 Mq 
(Smith ct al. 2007). The detection of soft unabsorbed X-rays by X MM — Newton indicated 
that the extended Circumstellar Matter (CSM) environment of SN 2006gy is of low density. 
Angoletto et al. (2008) challenged this hypothesis by considering highly opaque clumps 
distributed around SN 2006gy from which they derive a nickel mass estimate of about 3 Mq 
(see Smith et al. 2009 for an extensive summary). SN 2006gy might be consistent with a 
model of a pair formation supernova (Rakavy & Shaviv 1966; Barkat, Rakavy & Sack 1967; 
Smith et al. 2007), but SN 2005ap, which is even brighter but with a narrower light curve, 
cannot be (Quimby et al. 2007a). 

Recently, another apparently ultraluminous transient event was presented by the Su- 
pernovae Cosmology Project, SCP06F6 (Barbary et al. 2008). The possible high brightness, 
slow rise to maximum (100 days), and strange spectral features of SCP06F6 arc still under 
debate. There is no detected host galaxy consistent with the position of SCP06F6 although 
examples of low mass star-forming galaxies have been found that might be consistent with 
the upper limits (Dolphin et al. 2001). Thus SCP06F6 could define a new class if it is 
actually associated with a supernova explosion. Due to the peculiar spectral appearance of 
SCP06F6 and the lack of an identified host galaxy, its redshift remains uncertain. 

We discuss the uncertainties in the redshift, spectrum and luminosity in §2. In the 



-3- 



present work we examine a range of possibilities for the redshift and two special cases. We 
present a basic radioactive decay diffusion model plus the possible effects of a recombination 
front and dust formation in §3 and use it to place constraints on the distance, ejected mass, 
nickel mass, and kinetic energy in §4. In §5, we present models based on interaction of a 
supernova shock with circumstellar matter. A discussion is presented in §6. 

2. UNCERTAIN REDSHIFT, SPECTRUM AND LUMINOSITY 

Barbary et al. (2008) discussed the possibility of the extragalactic origin of SCP06F6 
since it has a small projected distance from the center of the galaxy cluster CL 1432.5+3332.8 
which has a redshift of 1.112. They noted that at this redshift, the absorption feature at 
5890 A aligns with Mgll AA 2796, 2803, but were unable to identify the other features. In 
particular, there are no features readily identifiable with hydrogen. In this work we argue 
that the redshift cannot be above 2; ~ 1 through various constraints on the ejecta mass, 
luminosity and radius. 

One special case corresponds to z = 0.143, as proposed by Gacnsikc ct al. (2008) 
by considering SCP06F6 as the explosion of a massive carbon-rich star with a cool and 
optically thick atmosphere. Gaensicke et al. (2008) propose that the broad "absorption" 
features present in the blue end of the spectrum of SCP06F6 (Figure 1) are molecular C2 
bands (the Swan bands) originating in the cool, optically-thick atmosphere of a carbon-rich 
progenitor. The presence of carbon molecules constrains the temperature of the atmosphere 
of the progenitor to be cool so that they are not destroyed. Gaensicke et al. (2008) find 
that carbon star spectra redshifted by 0.143 provide a decent fit to the combined VLT+Keck 
spectrum of SCP06F6 (see also Sokcr, Frankowski & Kashi 2008). Other features present 
in carbon star spectra, however, such as CN (A7900 A rest frame, Downes et 
al. 2004), do not seem to appeeir in the smoothed co-added spectra of SCP06F6 
(VLT-hSubaru-l-Keck). Other CN bands would fall too far to the blue to be 
constrained by the available spectral data. We return to this issue in §6. 

We also investigate the case for 2; ~ 0.57 where we can crudely identify three of the four 
broad absorption features by redshifting the spectrum of the local Type la SN 1992A or a 
template representing a Type IIP (Gililland et al. 1999) as seen in Figure 1. A redshift of 
z ~ 0.57 aligns the Ca II H&K absorption with the minimum at ~ 5900 A in SCP06F6. The 
two bluest minima crudely correspond to minima in the UV spectra of SN 1992A and the SN 
IIP template caused by overlapping lines of iron-peak elements (Harkness & Wheeler 1987; 
Kirshner et al. 1993). Similar features are seen in other supernovae of various spectral types 
with UV spectra (Bufano et al. 2009). The agreement is not good with the minimum at 
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about 5300 A in SCP06F6 where the SN II template has a peak. Note that at this redshift 
the continuum slopes redward of 4000 A (but not the spectral features) also roughly agree. 
The comparison in Figure 1 does not suggest that SCP06F6 is either a Type la or a Type II, 
but that SCP06F6 might have Ca II and iron-peak absorption in the rest-frame UV moving 
with velocities typical of supernovac. For reference, we also indicate in Figure 1 where Ha 
and H/3 would fall for this redshift. Ha would fall off the red end of the spectra and H/3 would 
fall near the gap in the spectra that is contaminated by telluric lines. For this redshift the 
lack of evidence for H in the spectrum cannot be taken as firm evidence for lack of hydrogen 
in the ejecta. The same is not true for the redshift of 0.143 where hydrogen might be present 
in emission or absorption for a variety of models, including those not based on supernovae. 

To approximate the bolometric lightcurve of SCP06F6 we use the VLT spectrum at 
maximum (Barbary et al. 2008) scaled to the i- and z- band fluxes observed at the same 
epoch. We then approximate the SED of the object with a triangle peaked at 6000 A that 
extends from 3000 to 12,000 A . We calculate the integral to get the quasi-bolometric flux at 
maximum light, F^oi- This approximate method gives Fhoi,max = 3.09 x 10~^^ erg s~^cm~^ 
which is in good agreement with the estimated value of 2.5 x 10~^^ erg s~^cm~'^ given in 
Barbary et al. (2008). We then scale the flux to the other epochs using the photometry in 
Table 1 of Barbary et al. (2008) assuming that the shape of the SED is roughly constant. 
We use the following formula to obtain bolometric luminosities: 

Lboi = Fboi^nr'^, (1) 

where r is the luminosity distance of the object which is equal to r = (1 -|- z/2) x cz/Hq, 
where c is the speed of hght, z is the redshift of the object and Hq — 70 km s~^Mpc~^ is 
the Hubble constant. This expression for r corresponds to an empty universe and is a good 
approximation to the full A-CDM expression. 



3. A SIMPLE DIFFUSION MODEL FOR SCP06F6 

3.1. Radioactive diffusion model 

We adopt the Arnett (1982) radioactive decay diffusion model as generalized by Valenti 
et al. (2008; see also Soderberg et al. 2008) to allow for both nickel and subsequent cobalt 
decay. We also allow for gamma-ray leakage. This simple model assumes a diffusive medium 
with a uniform density proflle and a photosphere that expands linearly in time. In general, 
the solutions we flnd for SCP06F6 imply large nickel and ejecta masses. With these assump- 
tions, we adopt the following expression for the diffusive luminosity released as a result of 
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radioactive decay: 

L{t) = M^.e-^'iiem - eco) f 2ze^'-'^ydz + eco f 2ze^'-'y^+'^^dz]{l - e"^*"'), (2) 

Jo Jo 

where x = t/tm, tm is the rise time to maximum (Arnett 2008), y = tm/2tNi with 
tm —8.8 days, s — tm{tco — iNi) /2tcotNi with tco =111-3 days, M^i is the initial nickel mass, 
e^vj = 3.9 X 10^° erg s~^g~^ and eco — 6.8 x 10^ erg s~^g~^ are the energy generation rates 
due to Ni and Co decay. The factor (1 — e~^^ ^) accounts for the gamma-ray leakage, where 
large A means that practically all gamma rays are trapped. The gamma-ray optical depth 
of the ejecta is taken to be = n^pR = At^"^, assuming spherical uniform density ejecta 
with radius R = vt and the Ni/Co confined in the center. This yields A = {SK^M^j) / {4:7rv^) 
which is controlled by the gamma-ray opacity /t^. The scaling follows from homologous 
expansion which is one of the basic assumptions of the simple analytic models that we adopt 
here. In general we find the effect of gamma-ray leakage to be small for the large masses 
implied in our models. 

Thus we have three main fitting parameters that determine the nature of the model 
light curve: the initial ^^Ni mass that determines its peak, the rise time to maximum tm 
that determines its width and A that determines the amount of gamma ray trapping and 
corresponds to a gamma-ray opacity. It is characteristic that, within this class of model, the 
nickel mass comes in only as a scaling parameter for the amplitude of the luminosity output. 
The more distant the source is, the brighter it is, since its intrinsic (rest frame) luminosity 
is higher compared to the observed luminosity. We thus expect a higher initial ^^Ni mass at 
higher redshifts. 

Unlike the maximum brightness, the width of the light curve in the object rest frame 
decreases with respect to its observed value, since the rest frame rise time, tm,rf, decreases 
with redshift: 

*"^"1 + Z' 

where tm,ob is the observed rise time. Assuming that the photosphere expands with approx- 
imately constant speed, Vph, we adopt the ejecta mass as a function of the observed rise 
time given by Valenti et al. (2008): 
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^-^'^ = lo^^-'^aTiF' 

where P is an integration constant equal to about 13.8 (Arnett 1982; Valenti et al. 2008), 
K is the mean opacity, and we have used Vph — {10 /SEke/ M^jY^'^ , as appropriate for the 
outer edge of a homologously expanding sphere of constant density with ejecta mass M^j and 
kinetic energy Eke- The ejecta mass is set by the rise time to maximum that depends on 
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the redshift of the object. Thus we obtain a scahng of the ejecta mass with redshift within 
the context of this model for SCP06F6. The photosphcric velocity of supernova ejecta is 
estimated from the width of lines in their spectra. For Type II supernovae it typically 
varies within the range 4000 — 6000 km s~^, and for Type la supernovae it is typically 
10, 000 — 15, 000 km s^^. We note that the features of SCP06F6 are generally broader than 
those of the Type IIP template of Figure 1 and perhaps more represented by the features of 
the Type la. Barbary et al. (2008) give Vph ~ 12, 000 km s'^. We adopt Vph = 10, 000 km 
as a fiducial value and assume this to refer to the photospheric velocity in the rest frame 
of the supernova for the models presented in §4. With assumed values of opacity we 
attempt to fit the observed light curve and to obtain constraints on physical properties of 
SCP06F6 that will in turn scale as a function of the redshift. 



3.2. The effects of recombination 

The hght curve presented by Barbary et al (2008) is roughly symmetric about the 
peak. The fact that SCP06F6 shows a symmetric hght curve is not in contradiction with 
observed supernovae characteristics. There are cases such as the Type lb supernova SN2005bf 
(Tominaga et al. 2005) for which the second peak has a very symmetric light curve. As we 
will show below, our simple radioactive diffusion models have difficulty accounting for the 
low values of the points 100 days after peak and hence reproducing the symmetric light 
curve. There are several factors that might account for this discrepancy. The fundamental 
assumption of the model of §3.1 is that the ejecta are optically thick. That may not be true 
at the last observed epoch. 

An expression for the optical depth at maximum light can be written as: 

Using the expression for nMej from Equation 4 and the kinematic radius Rk,max = Vphtmax/{^+ 
z) gives Tmax — c/vphi independent of redshift. For a photospheric velocity of 10,000 km s~^, 
the optical depth would be about 30, thus justifying the assumption of large optical depth at 
maximum light. Using the black body radius gives a similar estimate of large optical depth, 
but one that varies with redshift. 

As the ejecta expand the optical depth decreases as r ~ t^^. For constant opacity the 
optical depth at epoch 11 of Barbary et al. (2008) (Table 1) is thus Th ~ Tmax{tmax / tiiY — 
30/3 ~ 10 for Tmax — 30, tmax — 100 d and tn — 170 d. If the opacity dropped by a factor 
of several after maximum due to recombination in some element or elements, the light curve 
might decline more rapidly than our models. 
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To incorporate the effects of recombination we use the model developed by Arnett 
& Fu (1989) to fit the light curve of SN1987A. This model considers the effects of the 
recombination of hydrogen or helium. After peak luminosity the recombination front recedes 
inwards and the optical opacity drops dramatically. Thus the thermal energy generated by 
the presumed shock diffuses out more quickly and the decline is faster than it would be 
without recombination. The light curve that results from this model depends on 
the following input pcirameters: the initial radius of the progenitor Rq, the nickel 
mass yield M^i, the mass of the ejecta Mgj, the kinetic and the thermal energy 
Ekin and Eth respectively, the ionization temperature Tjon and the optical opacity 
K. The values that we used to obtain the fits shown in Figure 2 Eire summairized 
in Table 2. Also given in Table 2 is the mean ejecta velocity, Vmean — -s/^Ekin/M^ 
implied by the model fits. The amount of nickel needed to power the event in 
the recombination model is less than it is if the recombination effects are not 
taken into account, as expected. Specifically MNi,rec/Mjsfi ~ 0.6 for both redshift 
models, where Mjsn^rec is the mass of nickel that is required to power the peak of 
the light curve if the recombination effects are taken into account. 

As seen in Figure 2, a recombination model can produce a roughly symmetric light 
curve around the peak that subsequently follows the cobalt decay. In our model for 
SCP06F6 we consider the composition of the ejecta unknown and we adopt the 
Thompson scattering opacity values of 0.4 cm? at z =0.143, (a typical value 
for a H-rich envelope) and 0.2 cm^ at z =0.57 (a typical value for a H-poor 
envelope). As noted in §2 there is no evidence for hydrogen or helium in the spectrum of 
SCP06F6, but if it is at a redshift of 0.57, then the Ha line would be rcdshifted out of the 
range of the spectral converage by VLT and Subaru as shown in Barbary et al. (2008) and 
H/? would be contaminated by telluric lines (Figure 1). We are not able to accurately 
reproduce the light curve of SCP06F6, even when the recombination effects are 
taken into account, especially due to the fact that we cannot fit Epoch 11 of 
Bcirbary et al. (2008). It should be noted that the purpose of the recombination 
models that we consider is to illustrate the effects of this process on the light 
curve of a SN, and to show that recombination does tend to mcike the light 
curve more symmetric around the peak. More detailed models accounting for 
the density, composition and optical depth structure and appropriate radiative 
transfer cire beyond the scope of this work. 
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3.3. The effect of dust 

There could also be dust formation at late epochs. At Epoch 11 (Table 1, Barbary et 
al 2008) the observed flux is ~ 7 times lower than the flux that is predicted by the model 
(insensitive to redshift). Thus we can estimate the optical depth due to a hypothetical 
dust absorption using the formula Fn — Fo^nC"'^'*'"'*, where Fu is the flux at Epoch 11, 
Fo,ii is the flux at epoch 11 as predicted by the radioactive decay diffusion model and Tdust 
is the dust optical depth. This yields Tdust — 1-9 if we consider a spherical, optically thin 
dust shell at the outer edge of the SN ejecta with radius ~ Rn = Vphtu and thickness 
dR « Rii. The required dust mass is then Mdust,ii — ^T^R\iPdustdR, and the dust optical 
depth Tdust = tidustPdustdR. Combining these equations yields Mdust = r / {Kdust'^T<:Rli) where 
f^dust is the opacity of the dust grains. We adopt the value Kdust = 10^ cm^ g^^ as a 
gross estimate, using the Draine (2003) model extinction curves for Milky way dust grain 
composition. For this choice of parameters we estimate the dust mass to be ~ 3 x 10~^ Mq. 
This value is in agreement with values that are estimated for other Type-II supernovae 
(Kotak et al. 2005; Kotak et al. 2009). 

If the dust formation took place well inside the ejecta, for example at u ~ 1000 km s^^, 
as in Type II-P SNe (Kotak et al. 2009), then the dust mass would be 100 times higher than 
the value estimated above, on the order of 0.02Mq. This is somewhat higher than has been 
observed in other SNe, but still not unrealistic, taking into account that the dust masses 
derived in this way are only lower limits. If the dust distribution is not uniform and optically 
thin but rather it is distributed in optically thick clumps (Ercolano et al., 2007), then the 
total dust mass could be much higher than that derived from the optically-thin model. 

Note that the energy absorbed by dust is expected to be re-radiated as mid-IR photons 
after the dust particles (or clumps) are heated up to Tdust ~ 600 - 1000 K. Thus, strictly 
speaking, the true bolometric luminosity remains unchanged in this case. In §2, we estimated 
the bolometric luminosity by assuming that all radiation comes from the optical and near-IR 
regime. The Arnett model (Equation 2) also uses essentially the same assumption due to 
the thermalization of the gamma-rays. Within this context the optical photons absorbed 
by dust can be considered as "lost," even though their energy is expected to re-appear as a 
mid-IR excess, outside the original SED. 
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4. CONSTRAINTS ON SCP06F6 

4.1. General redshift constraints 

Here we calculate the redshift dependence of the physical petrameters of the 
Arnett model given in Equation 2, without recombination. Figure 3 shows how 
the ejecta mass and the nickel mass depend on redshift based on fits to the light curve 
for various choices of the mean opacity. The recombination models have the same 
redshift dependence, but the nickel masses are about 40% lower, the ejecta 
masses are about 30% lower and the kinetic energies are lower by a factor of 2 to 
3 compared to our fiducial models with the same opacity. The ejecta mass for the 
basic radioactive diff"usion models is calculated using Equation 4 adopting tm,ob = 100 d 
for the rise time to maximum in the observer's frame and adopting Vph = 10, 000 km 
as the fiducial velocity. The nickel mass is calculated by assuming tmax — tm and solving 
Equation 2 for M^i given the observed value of maximum luminosity, which scales with 
redshift according to Equation 2. We consider values of the mean opacity ranging from 
0.05 cm^ g'^ as might be representative of a metal-rich composition (Sutherland & Wheeler 
1984; Soderberg et al. 2008) up to 0.40 crn^ as might be representative of electron 
scattering in a pure hydrogen plasma. Figure 3 also shows the best fitting results for the 
nickel mass and the ejecta mass in the case of k=0.1 an^ g^^ for the same redshift range 
and with redshift increments of 0.1 (filled circles and squares respectively). To determine 
the best fitting parameters {M^i, tm and A) for each redshift we developed a simple Monte 
Carlo chi-square minimization code that scans through all the parameter space and finds 
the minimum value and the parameters that correspond to that value. As can be seen in 
Figure 3, the best fitting results are in very good agreement with the analytic ones assuming 
tm,ob — 100 d as a fixed petrameter. 

The region above redshift of ~1.1 is forbidden since the required nickel mass would 
exceed the required ejecta mass to account for the width of the light curve even for rather 
small opacity. The region to the right of the dashed vertical line at z ~0.8 (for Ae=0.1 cm~^ 
g~^) is also forbidden in practice because it is unlikely that the nickel mass exceeds half 
the total ejecta mass in an astrophysically realistic model of, for example, a pair-formation 
or core-collapse supernova. These redshift limits get tighter for larger mean opacities or 
lower photospheric velocities. For electron scattering in a pure ionized hydrogen plasma, the 
redshift would have to be less than about 0.65 for this sort of model to be self-consistent. In- 
voking recombination does not change this constraint substantially since Mej/Mm 
is about the same (Tables 1,2). A nickel mass of IMq that would be characteristic of 
a Type la is obtained at a redshift of about 0.3. At this redshift, the ejecta would exceed 
30 Mq for a mean opacity of k ^ 0.1 cm^ g~^, so the explosion is certainly not that of an 
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exploding white dwarf. 

Figure 4 shows the scahng with redshift of the kinetic energy, E^e = V^-^ej^mean' 
and the total radiated energy, E^ad- The kinetic energy is determined from the ejecta mass, 
assuming a mean velocity of Vmean — ^/^/5vph — 7, 800 km and the same range of mean 
opacity as before. Note that the kinetic energy is especially sensitive to the assumed photo- 
spheric velocity, scaling as Eke oc (from Equation 4). The more tightly-constrained 
recombination models require a lower Ekin/M^j and hence Vmean than our fiducial 
radioactive diffusion models (Table 2). This raises the possibility that our fidu- 
cial models presented here and in Figure 4 overestimate the ejecta velocity and 
kinetic energy. The radiated energy is obtained by integrating the bolometric luminosity 
over the whole light curve. Since a strong shock distributes energy equally between kinetic 
energy and thermal energy in its wake and some thermal energy is likely to be subsequently 
converted to kinetic energy by adiabatic expansion and PdV work, it is unlikely that the 
radiated energy can exceed the kinetic energy. This constraint gives another limit on the 
redshift, z ^ 1.1 for k.=0.1 cm~^ g~^, as shown in Figure 3. Note that if the opacity were 
larger or the photospheric velocity smaller, this constraint also gets tighter. 

An estimate of the radius of the photosphere can be obtained from kinematics {R — Vpht) 
or emissivity (a blackbody radius). Figure 5 gives the radius estimated at maximum light in 
these two ways assuming a constant rest frame photospheric velocity of 10,000 km s^^ and 
a constant rest frame blackbody temperature of 5000 K. Note that these two estimates only 
agree at a single redshift. Each estimate is uncertain, so we also denote the range over which 
the radii agree within a factor of two or a factor of three. The latter restricts the range of 
redshift to z ~ 0.2 - 0.9. 

Another consistency check follows from noticing that the slope of the continuum varies 
little over the three epochs for which spectra were given by Barbary ct al. (2008). We 
have estimated the effective temperature in our models at these epochs using the bolometric 
luminosity and the kinematic radius from Figure 5. The effective temperature derived this 
way (that is formally 5000 K at z~0.5) varies by ^ 300 K over the spectral epochs, basically 
consistent with the observations. 



4.2. Special case 1: z — 0.143 

We next examine the special case of z = 0.143 as suggested by Gaensicke et al. (2008) 
who fit the spectrum of SCP06F6 with redshifted carbon-rich stellar spectra from the SDSS 
database. The result of the analysis for z = 0.143 is shown in the top panel of Figure 2. The 
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best fit gives a value of Mm = 0.28 Mq and Me,=37.74 ( q.i cm%-i )~' ( io,ooo\ms-i ) ^0- 
The best fitting A parameter which controls the gamma ray leakage implies a small gamma 
ray opacity = 5 x 10~^ cm^ g~^, whereas a typical gamma-ray opacity might be ~0.03 
cm^ ( Colgate, Petschek and Kriese 1980). Gamma-ray leakage is not an important effect 
for the large masses derived. The initial nickel mass is within the range expected for typical 
SNe while the ejecta mass is rather large, due to the fact that SCP06F6 shows a significantly 
slow rise to maximum and a slow dechne. Note that the ejecta mass varies inversely with the 
assumed, and uncertain, opacity (Equation 4). This mass could be less if the mean opacity 
were greater than the value, k, — 0.1 cm"^ g~^ assumed here for illustration. It would also 
be less for smaller Vphot- 

The value of the ejecta mass implies a total kinetic energy of Eke = 2.4 x 10^^ 

( o.i cm-^ g-i )~^ ( io,ooo\m s-i )^ ^ radiated energy is 9.1 x 10^^ 

erg. While the estimated radiated energy is roughly in accord with "normal" supernovae, 
the estimated kinetic energy is very large. For the assumed opacity, the implication is that 
the initial radius was rather small compared to the radius of the photosphere at maximum 
light so that substantial initial thermal energy was lost to adiabatic expansion. For a larger 
opacity and a photospheric velocity more in the range typical of Type 11 supernovae, the 
estimated kinetic and radiated energies would be more nearly equivalent. 

The top panel of Figure 2 also shows the results of a model that incorporates the 
effects of recombination on the light curve (§3.2) for this redshift case. The recombination 
temperature was found to be T = 5500 K with k =0.4 cm^ g^^ which is appropriate 
for a H-rich atmosphere, although there is no evidence for the presence of H in 
the spectrum if this were the correct redshift. For this model the nickel mass, 
the ejecta mass and the kinetic energy cire lower than for the basic model, as 
mentioned above. The kinetic energy, and hence implicitly the mean velocity, was varied 
to produce the dashed curve in the top panel of Figure 2. Given the uncertainties, no 
attempt was made to determine the "best fit" over the full parameter range. This model 
shows that recombination could help to produce a roughly symmetric light curve and that 
this physics, in principle, is relevant to SCP06F6, although it was not possible to fit 
the last measured data point (Epoch 11 of Barbary et al.). We note that for this 
redshift, the lack of H or He features in the spectrum is an issue (§2). We also note that 
dust formation may play a role in the decline from peak (§3.3). 
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4.3. Special case 2: z — 0.57 

We now follow the same argument as before for the case of redshift of 0.57. The results 
of this fit are shown in the bottom panel of Figure 2. The best fit gives an estimate of 
M^i = 4.83 Mq and Me, = 19.68 ( q.i cm-2g-i )~' ( io,ooo\m s-O ^0- ^ote that the ejecta mass 
is less in this case than for z=0.143 because the dilation reduces the width of the 
Ught curve in the rest frame. Again the ejecta mass is rather large, but scales with the 
uncertain opacity and velocity. The nickel mass in this case is higher, as expected, but within 
the range of values predicted for other luminous SNe based on radioactive diffusion models 
(Quimby et al. 2007; Smith et al. 2007). The best fitting gamma-ray leakage parameter in 
this case yields a gamma ray opacity of k,^ = 0.03 cm^ which is in agreement with the 
generally assumed value of 0.03 cm^ g~^. Gamma-ray leakage has no substantial effect on 
the light curve. 

The total ejecta mass in this case implies a kinetic energy equal to Eke = 1-2 x 10^^ 
( o.i cm-2 g-i )~^ ( lo.ooo'km s'l )^ ^^S, with a total radiated energy of Erad = 1-4 x 10^^ erg. 
While the radiated energy is somewhat large for a normal supernova, the kinetic energy is 
again quite large, and comparable to that estimated for z = 0.143. For the fiducial parameters 
in this case, the implication is again that much of the initial shock energy must have been 
lost to adiabatic expansion. If the opacity were larger than adopted here for illustration and 
the photospheric velocity somewhat smaller, as impUed by the recombination models, 
the estimated kinetic energy could be more representative of normal supernovae and the 
estimated radiated energy could be comparable to the estimated kinetic energy, in which 
case substantial adiabatic losses would not be implied. This would imply a large initial 
radius and thus perhaps a dense circumstellar medium as will be discussed below. 

As for the lower redshift case, the basic radioactive diffusion model does not produce 
the steep post-maximum decline. The dashed line in the bottom panel of Figure 2 shows a 
model with recombination at 10,000 K and k =0.2 cm^ which is appropriate for a 
H-poor envelope (see Table 2). Dust formation might also play a role (§3.3). We expect 
no evidence for H to be observed at this redshift (§2). 

5. CSM INTERACTION 

Now we investigate the possibility of the contribution of CSM interaction producing 
optical/NIR emission in the observed light curve of SCP06F6. There are two versions of 
this circumstance. In the model of Smith & McCray (2007), the CSM shell is optically thick 
and the light curve is controlled by diffusion. Alternatively, the CSM could be optically thin 
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enough to radiate the shock coUision energy "instantaneously," but optically thick enough 
to convert the shock energy to optical radiation. 

In the case where the energy release is dominated by diffusion, the energy to power 
the light curve arises in the putative collision of the supernova ejecta with the extended 
dense CSM. There is no need for any radioactive input to power the light curve 
near maximum light, although some contribution from that source cannot be 
ruled out. This class of model gives no simple way to estimate the maximum luminosity in 
terms of physical input parameters. Those input parameters would be the initial radius and 
kinetic energy of the underlying supernova and the initial radius and density distribution 
of the CSM. An extreme version of this class of model is obtained with the assumption 
that the initial radius of the circumstellar shell was not far from the observed radius of 
the photosphere at maximum light so that adiabatic losses within the shell are assumed 
to be minimal. In this case, the post-shock thermal energy content of the CSM envelope 
should be roughly comparable to the radiated energy and to the kinetic energy and both 
should be roughly comparable to the initial supernova shock energy, the kinetic energy of 
the underlying supernova (within factors of two). These aspects mean that, except for this 
"extreme" version, there is no simple way to determine a physical parameter analogous to the 
initial nickel mass in the radioactive diffusion model derived from the light curve peak. On 
the other hand, this model still assumes that the luminous output derives from the diffusion 
of thermal energy from within the optically-thick ejecta. In this case, the mass of the ejecta 
that determines the width of the light curve is roughly the sum of the underlying supernova 
ejecta and the mass of the shocked CSM matter (Smith & McCray 2007), and the tools 
presented in §3 still allow a determination of this effective mass as a function of redshift. 

If the optical depth is small enough that the diffusion time for the release of post-shock 
energy is short, then the luminosity released as the CSM is shocked is given by (Ofek et al. 
2007; Gezari et al. 2008; Smith et al. 2008): 

L ~ 2T:pcsMR^vlf, (6) 

where pcsM is the density of the CSM and Vgh is the velocity of the shock. Note 
that in this case, the luminosity is presumed to reflect the local density, and hence the shape 
of the light curve is given by the density distribution. The light curve can be reproduced by 
a suitable, although entirely ad hoc, assumption of a density profile. We also note that while 
some information on the mass ejected into the CSM can be estimated, this model provides 
no useful separate information on the ejecta mass analogous to the constraint of the rise 
time to maximum. There is thus also no independent constraint on the associated kinetic 
energy. 
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For the CSM diffusion models, the width of the hght curve yields an estimate of the 
diffusion time and hence the ejecta mass from Equation 4. An estimate of the radius of the 
shell from kinematics or emissivity (Figure 5) then yields an estimate of the mean density, 
p. For the z=0.143 case the radius at maximum light is estimated from the kinematics to 
be -Rfc.0.143 = Vphtrnax/i)- + z) = 7.5 X 10^^ cm for a photospheric velocity of 10,000 km . 
Assuming a temperature of 5000 K in the rest frame of the supernova (Gaensicke et al., 
2008), the black body radius at maximum light is Rbb,o.i43, = 1-8 x 10^^ cm. The rather large 
discrepancy in these two methods for estimating the radius near maximum light suggests 
a basic inconsistency in the diffusion models for this smaller redshift. Following the same 
method for the z=0.57 case, we find that the blackbody radius is Rbb,o.57 — 6.2 x 10^^ cm while 
the kinematic radius at maximum light is Rk.0.57 = 5.5 x 10^^ cm. This represents generally 
good agreement. If the true blackbody radius at maximum light were much higher, this 
model might be strained. 

The significantly low fraction of radiation energy compared to the kinetic energy for the 
model at z = 0.143, Eth/Ekin ~ 0.005 (§4.3), implies that there must have been substantial 
adiabatic losses for this CSM diffusion model to be viable, unless the opacity is substantially 
larger than the fiducial value of k =0.1 cm^ g^^ (sec Figure 4) and the photospheric velocity 
substantially lower. For the case at z = 0.57, the radiated energy is again estimated to be 
small compared to the kinetic energy, but rather modest changes in the fiducial parameters 
(larger opacity, smaller velocity) could make the radiated energy comparable to the kinetic 
energy. In this case, the light curve could be accounted for by the collision of an underlying 
supernova with a large, dense, CSM shell. For this situation, with negligible adiabatic losses, 
the energy of the underlying supernova would be about twice the radiated energy, or about 
3 x 10^^ ergs. The estimate of optical depth near maximum light, Tmax — 30, from Equation 
5 also applies to this diffusion model. This constant value of opacity is given as the horizontal 
dashed line in Figure 6. Because a CSM envelope is likely to be of relatively large opacity, 
the shock diffusion model may be appropriate to this redshift. The expected composition 
of the CSM shell is hydrogen so the lack for H features in the spectrum is a constraint at 
z=0.143 but not necessarily at z=0.57 (§2). 

An estimate of the optical depth can also be made in the context of the model for which 
the optical depth is modest and the shock energy is radiated "instantaneously." In this case, 
we can estimate the density with which the shock collides near maximum light from Equation 
5 as: ^ 

^" ^max^sh 

Taking the shock velocity to be constant and equal to the photospheric expansion velocity 
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and hence Umax = Vphtmax yields: 



'max 



(8) 




with tmax — tm,obs/{^ + z) . Note that the optical depth for this model is very sensitive 
to the assumed shock velocity. Adopting a characteristic shock velocity of 10,000 km s~^, 
we can estimate r for a given opacity. In Figure 6 we also show how the optical depth 
at maximum light estimated in this way scales with redshift. Note that the optical depth 
estimated based on this model increases rapidly with redshift. For optical depths less than 
unity, the luminosity is expected to be emitted at high-energy, non-optical wavelengths. For 
optical depths much greater that unity, the diffusion time, not shock propagation, will control 
the emission timescale. From Figure 6, the range of validity of this model for the fiducial 
parameters is z > 0.9. For the z = 0.57 case, this model could be made self-consistent for 
higher opacity or smaller shock velocities. The lack of evidence for H in the spectrum again 
constrains this class of models, especially at lower redshifts. 

We can also constrain the mass loss rate associated with the CSM. For circumstances 
in which the CSM is optically-thick and diffusion controls the light curve, we can use the 
estimates of the CSM mass and the radius of the configuration to make an estimate of the 
effective mass loss rate. The mass loss rate can be estimated as: 



where v^, is the velocity of the "wind" that led to the formation of the CSM. For a red-giant 
type wind, this velocity might be 10 km s~^. For an LBV- type mass loss event, a typical 
velocity might be ~ 100 km s"^ (Smith et al. 2004; Smith 2006). If we assume that the 
diffusion time is dominated by the mass ejected into the CSM prior to explosion, and hence 
that the "ejected" mass in Figure 2 is a measure of the CSM mass, then we can estimate the 
effective mass loss rate. For z = 0.143, we get 0.15 M© yr^^ and 1.5 Mq yr~^ for equal 
to 10 and 100 km s~\ respectively. For z = 0.57, we get 0.11 M0 yr~^ and 1.11 Mq yr~^. 
These values were all based on our fiducial opacity and so scale as ( q 1 cm-2 g-i )""*"• ^ 
Wolf-Rayet progenitor the wind velocity could be up to 1000 km s~^ implying an even larger 
mass loss rate. All these values are large, implying a more LBV-fike process if this diffusion 
model is pertinent. 

We can also estimate the mass loss rate implied in the model where the luminosity arises 
by shock interaction in a CSM of modest optical depth. The total envelope mass within a 
radius R — Vshtmax is, with Equation 7: 




(9) 



2L 




(10) 
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With Equation 9, we can then write: 

M=-— 3— -10 MqVT —3 , (11) 

^ "sh '^s?i,10,000 

where I/maa;,42 IS the peak luminosity in 10^^ erg s~^, t'u,,ioo is the wind velocity in units of 100 
km and Vsh,io,ooo is the shock speed in units of 10,000 km s~^. For a redshift of 0.57, the 
corresponding mass loss rate is about 3 x 10~^Vyj^ioo '^sh loooo yr~^- This mass loss rate 
could be representative of a Wolf-Rayet star progenitor, but then the wind velocity should 
be higher to correspond to observed Wolf-Rayet stars. Another difficulty is that this mass 
distribution cannot be a standard p cx wind as roughly expected for a Wolf-Rayet star. 
Rather, the density profile must be carefully "designed" to reproduce the shape of the light 
curve. This conceptual problem mitigates against this model on general grounds. 

Gaensicke et al. (2008) reported a very high 0.5 - 10 keV X-ray flux near optical 
maximum of about 10~^^ erg s^^ cm~^. This is nearly foiu times the bolomctric flux. It is 
possible that the reported detection should instead be treated as an upper limit (D. Pooley, 
private communication). The X-ray flux could arise from shocked optically-thin CSM. In 
this case Equations 5 and 7 pertain, and if we take the reported X-ray flux at face value, we 
estimate a particle density 6.4 x 10^ cm~^ and 2.6 x 10^ cm~^ for redshifts of 0.143 and 0.57, 
respectively, for a shock velocity of 10,000 km s^^. This density is 100-10,000 times smaller 
than the average density p of a dense optically-thick shell for the two redshift cases. This 
result implies that the observed X-ray flux, if real, might possibly arise in the outer, lower 
density portion of the diffusive shell that produces the optical display. For a CSM of modest 
optical depth, it is difficult to see how a self-consistent model based on shock interaction 
could simultaneously explain the optical and the X-ray luminosity. A possibility is that the 
medium is clumpy with dense clumps providing the optical emission and inter-clump regions 
providing X-rays. 



6. DISCUSSION AND SUMMARY 

We discussed the applicability of some widely used SN light curve models in reproducing 
the observed light curve of the luminous peculiar transient SCP06F6 discovered by Barbary 
et al. (2008). The parameters estimated based on various models arc summarized in Table 1. 
The observed light curve can be approximated by a smooth simple diffusion model, but the 
rapid decline from maximum is not consistent with the simple models. Fits to a radioactive 
decay diffusion model provide estimates for the ejecta mass, the nickel mass, and the gamma 
ray opacity versus redshift depending on key, but uncertain parameters, especially the mean 
optical opacity and photospheric velocity. We also note that while the data invite the 
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interpretation of a smooth rise and decline of the hght curve, the photometric data are sparse. 
The possibihty of a precursor peak, such as displayed by SN1987A or SN2005bf, cannot be 
ruled out. The existence of such a feature would modify any interpretation of the data. 

We considered two specific choices of the redshift, z — 0.143 based on a suggestion 
by Gaensicke et al. (2008) and z — 0.57, based on a crude fit to Ca H&K and iron peak 
absorptions in various supernovae. Substantially higher rcdshifts, greater than ~ 1.1, do 
not lead to self-consisent results in the context of the radioactive decay diffusion model: 
the nickel mass becomes too large compared to the ejecta mass and the radiated energy 
exceeds the kinetic energy. For these two specific choices of redshift, wc found nickel masses 
of 0.28 and 4.83 M©, respectively, and ejecta masses of ~37 and ~20 M©, respectively, for 
an adopted mean opacity of 0.1 err? g'^ and a photospheric velocity of 10,000 km s~^. The 
estimated ejecta mass scales inversely with the opacity and directly with the velocity. The 
results are roughly commensurate with normal core-collapse supernovae at lower redshift, 
but at the larger redshift, the nickel mass would require a different situation, perhaps similar 
to that invoked for some models of SN 2006gy based on pair-formation supernovae (Smith 
et al. 2007). For both redshifts, the implied kinetic energy is very high, ~ 10^^ erg s~^ for 
the fidicial parameters, but would be less for smaller ejecta velocities. In either case, 
our simple diffusion models do not acount for the rapid post maximum decline. Some other 
factor would need to be invoked such as a change in opacity or dust formation. We show 
that recombination of H, maybe He, might plausibly account for a decline in opacity and 
hence in the light curve. 

We also consider models in which the optical luminosity is provided by collision of the 
supernova with a dense GSM shell in the spirit of Smith & McCray (2008) or models in which 
the shock energy is thermalized, but rapidly radiated away (Ofek et al. 2007; Gezari et al. 
2008; Smith et al. 2008). In general, the models with dense, optically thick extended shells 
behave in a manner similar to that of the radioactive decay diffusion models, except that 
the energy to heat the GSM envelope is presumed to derive from the kinetic energy of the 
underlying explosion, rather than radioactive decay. These models applied at low redshift 
give the radiated energy substantially less than the kinetic energy, implying that to be self- 
consistent, the GSM envelope must have undergone substantial adiabatic expansion after 
being shocked. The same is true for the model apphed at larger redshifts, but reasonable 
choices of the parameters would yield a model for which the shell had expanded httle by 
maximum light. For such choice of parameters, the energy of the underlying explosion could 
be more modest. These models demand substantial mass in the GSM and hence effective 
mass loss rates that are reminiscent of LBV mass-loss episodes. Models in which the GSM 
is only modestly optically thick so the energy can be thermalized but also radiated rapidly 
could be applicable at higher redshifts. Such models might be consistent with mass loss 
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rates reminiscent of Wolf-Rayet stars, but to match the shape of the hght curve, the density 
profile would need to be carefully contrived, an unlikely happenstance. 

The spectrum of SCP06F6 does not resemble any other transient spectrum that has 
been obtained so far. The four broad absorption features can be reasonably well fit by the 
molecular C2 Swan bands of a carbon-rich star of temperature 3000-10, OOOK (Gaensicke et al. 
2008). Swan bands are also observed in white dwarfs at substantially higher temperatures, 
but only because of the very high gravity that would not pertain to supernova ejecta. Pending 
the calculation of a realistic supernova atmosphere model showing that C2 can form and 
survive, we find this hypothesis intriguing, but unlikely, at least around maximum light 
when the ejecta is supposted to have Tg// ~ 5000K. On the other hand, at later phases, 
when the ejecta has cooled, the formation of C2 and other molecules, even dust grains, is 
indeed a possibihty although this would imply the presence of CN (A7900 A rest 
frame) which we do not detect in the co-added VLT, Subeiru and Keck spectra 
of SCP06F6. 

We propose that a model of SCP06F6 worth further consideration is one in which the 
event is at a redshift of z ~ 0.57 for which the blue absorption features are Ca H&K and 
iron-peak absorption features as qualitatively seen in the UV of some nearby supernovae. At 
this redshift the light curve can be reproduced semi- quantitatively with a diffusion model 
based either on radioactive decay requiring several solar masses of nickel or collision of a 
moderately energetic supernova with a dense GSM envelope. 

The characteristics of SGP06F6 that are deduced from its hght curve at z—0.57 indicate 
that its luminosity is one order of magnitude less than that of the exceptionally luminous 
2005ap, 2006gy, 2008es and 2005tf. This means that if SCP06F6 is a supernova fitting the 
characteristics presented here, defining a new class of objects is not necessary to account for 
it. It fits well within the rubric of these other bright events. 

We are grateful to Kyle Barbary, Robert Quimby, Andy Howell, Nick Suntzeff, Alicia 
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Fig. 1. — Comparison of the VLT and the Subaru observer's frame spectrum of SCP 06F6 
obtained on 05-18-2006 (near maximum) and 05-22-2006 respectively (Barbary et al. 2008) 
with the lUE+CTIO spectrum of the Type la SN 1992A obtained on 01-24-1992 (+18d after 
maximum) (Kirshner et ah 1993) and with a template Type-IIP spectrum at -|-6d after the 
explosion (Gilliland et al. 1999) both boosted to a redshift of 0.57. The dashed vertical hnes 
indicate the positions of the Ca H&K absorption component and the Mg II triplet. The 
ranges of Fe and Si blends are also indicated. Line identification is based on Kirshner et al. 
(1993). The redshifted positions of the Ha and H/9 lines are also indicated for illustration. 
Note that Ha and H/3 could be present, but unobserved, at this redshift due to the large 
redshift and telluric contamination. 
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Fig. 2. — The rest frame quasi-bolometric luminosity light curve of SCP 06F6 (§3.1) at 
redshifts of 0.143 (top panel) and 0.57 (bottom panel) (sohd points). In each case, the solid 
line is a simple radioactive diffusion model (§3.1) and the dashed line an illustration of the 
Arnett & Fu (1989) model that includes the effects of recombination giving a decline in the 
optical opacity and thus a rapid post-maximum decline. The derived values of the original 
nickel mass and of the ejecta mass are given in the insets. Note that in this figure 
the basic model assumes k = 0.1 cm^ and the corresponding recombination 
models assume k = 0.4 and 0.2 cm^ g \ respectively, for z — 0.143 and z = 0.57 
(see Tables 1, 2). 
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Fig. 3. — The dependence of the initial nickel mass (rising line) and the ejecta mass (decreas- 
ing lines) on redshift for SCP06F6. The line for the nickel mass is based on Equations 1 and 
2. The lines for the ejecta mass are from Equation 4. The ejecta mass depends linearly with 
the photospheric velocity assumed to be 10,000 km s~^ and inversely on the mean opacity, 
for which three choices are shown: n = 0.05 (solid) ,0.1 (dotted dashed) ,0.4 (dotted) cm^ 
g~^. The filled circles and squares correspond to the best fitting estimates of the nickel 
mass and the ejecta mass respectively calculated by the Monte-Carlo minimization code 
(see text). At redshifts to the right of the solid vertical line the mass of radioactive nickel 
becomes larger than the total ejecta mass for k = 0.1 cm^ g~^. This line thus defines the 
lower redshift boundary of a "forbidden region" in redshift space for this class of radioactive 
diffusion model. The dashed vertical line indicates the redshift at which the nickel mass 
is half the total ejected mass for k = 0.1 cm^ g^^. The forbidden region extends to lower 
redshift for higher opacity and lower photospheric velocity. 
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Fig. 4. — The dependence of the total radiated energy (thick sohd hne) and the kinetic 
energy of the ejecta on redshift for SCP06F6. Three choices of the mean opacity are shown 
for the kinetic energy (see Figure 1). The region at redshifts to the right of the dashed 
vertical line at 2; = 1.18 is forbidden for k, — 0.1 cm^ g~^ since the radiated energy becomes 
larger than the kinetic energy of the ejecta. This forbidden region extends to lower redshifts 
for larger opacity. 
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Fig. 5. — The dependence on redshift of the photospheric radii based on kinematics (R = Yph 
t) and on black body emission at maximum hght. The kinematic radius assumes a constant 
photospheric velocity in the rest frame of 10,000 km s~^. The black body radius assumes a 
constant temperature of 5000 K in the rest frame. The vertical lines indicate the redshift 
range over which these estimates of the radius agree within a factor of two (solid line) and 
a factor of three (dashed hne) , respectively. 
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Fig. 6. — The dependence of the optical depth at maximum light on redshift for CSM models 
of SCP06F6. The horizontal dot-dash line at r = 30 represents the estimated optical depth 
(r = c/vph) for a diffusion shell model with radius given by the kinematic estimate (see text). 
The solid horizontal line corresponds to r = 1. The thin solid, dashed and dotted lines are 
based on a model of moderate optical depth for which the luminosity is quickly radiated 
behind the shock. This model is valid only for t^I. 
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Table 1. Summary of the fiducial physical parameters of SCP06F6 







z = 0.143 






z = 0.57 




k/ erring ^ 


0.05 


0.1 


0.4 


0.05 


0.1 


0.4 


*m,7-//rfoi/s 


90 


90 


90 


65 


65 


65 




0.28 


0.28 


0.28 


4.83 


4.83 


4.83 




75.4 


37.7 


9.4 


39.4 


19.7 


4.9 


Bth/lO^Oergs 


0.91 


0.91 


0.91 


14.0 


14.0 


14.0 




42.0 


21.0 


5.3 


5.5 


11.0 


2.8 


fli,i,/10l^cm 


1.7 


1.7 


1.7 


7.8 


7.8 


7.8 




7.5 


7.5 


7.5 


5.5 


5.5 


5.5 




30 


30 


30 


30 


30 


30 




0.001 


0.003 


0.011 


0.04 


0.08 


0.32 


Mthick/MiT, s/r-lf 


3.0 


1.5 


0.38 


2.22 


1.11 


0.28 




2 X IQ-" 


2 X IQ-" 


2 X 10-" 


3 X 10-^ 


3 X 10-^ 


3 X 10-^ 



Note. — (vp/i = 10,000 kms ^ for all cases). fThe values of Tf^ick ^thick estimated based 

on a model of an optically thick shell and the values of Tffj_j^^ and M^f^^^ are estimated based on a model 
of moderate optical depth for which the luminosity is quickly radiated behind the shock. The given mass 
loss rates assume a "wind" velocity of 100 km s""*^. 



Table 2. Summary of the parameters of the recombination models 



Parameter 


z = 0.143 


z = 0.57 


iJo/lO^^cm 


2 


2 




0.16 


3.1 




6.2 


6.9 


Bth/loSlergs 


1.5 


1.0 




1.5 


3.0 


Tion./K 


5500 


10000 


K /cm^ff — ^ 


0.4 


0.2 




3,900 


2,700 



